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QCD phase diagram is obtained by analysing centrality dependence of J/t/) suppression in 
y^^^ = 17.3 Pb+Pb and ^/sj^j^— 200 GeV Au+Au collisions. J/ip's produced in initial interac- 
tions are assumed to dissolve if local temperature exceeds a threshold temperature. The threshold 
temperature depends on the (local) fluid temperature and baryonic chemical potential, which are 
obtained, under certain assumptions, from experimentally determined quantities e.g. rapidity den- 
sity, net baryon density. QCD critical line with curvature parameter k ~ 0.018 — 0.23 is consistent 
with experimental J/tp data, k is factor of ~3 larger than in lattice QCD calculations. 
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I. INTRODUCTION 

In recent years, there is much interest in QGP phase 
diagram in T — /i plane, T the temperature and fj, the 
baryonic chemical potential [H,[l])lll)lll- Presently, only 
one point in the phase diagram is known from lattice 
QCD simulations, at /i «0, the QCD transition is a cross 
over ^5|| at (pseudo) critical temperature T = wlS? 
(3) (3) MeV p. From theoretical considerations, QCD 
phase transition is expected to be 1st order in baryon 
dense matter. One then expects a QCD critical end point 
(CEP), where the transition changes from 1st order to 
cross-over. To access the QCD phase diagram and CEP, 
systematic study of nuclear collisions, as a function of col- 
lision energy, has been planned at FAIR Q in the energy 
range i?zab=5-40 GeV. 

At finite baryon density, Fermion determinant is com- 
plex and standard technique of Monte-Carlo importance 
sampling fails. Several techniques have been suggested to 
circumvent the problem, (i) reweighting 01,0 , (h) ana- 
lytical continuation of imaginary chemical potential [lo| , 
[llj and (iii) Taylor expansion 

g3]>[Il- These methods 
has been used to locate the phase boundary in T — /i 
plane. The calculations suggest that the curvature pa- 
rameter [k) in the expansion, 



is small As an example, in Figll] QCD phase dia- 
gram obtained in the analytical continuation method [8] 
(the filled circles) and in Taylor expansion pJi] of (chi- 
ral) are shown. Both the methods gives nearly identical 
phase diagram for ^/Tdp — 0) < 3GeV, curvature pa- 
rameter is small, k w 0.006. At larger /i, they differ 
marginally. For comparison, in Figll] chemical freeze- 
out curve [l^,[l3|; obtained in statistical model analysis 
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of particle ratios are shown (the red line). Curvature of 
the chemical freeze-out curve is factor of 4 larger than 
the curvature in the QCD phase diagram. Small curva- 
ture of the QCD critical line, compared to the chemi- 
cal freeze-out is interesting. Experimental signal of crit- 
ical end point will get diluted as the deconfined medium 
produced at the critical end point will evolve longer to 
reach chemical freeze-out. Fluid will have more time to 
washout any signature of CEP. 

In the present letter, we obtain the curvature parame- 
ter K by analysing experimental data on J/ip suppres- 
sion. Existing data in y^j^^=17.3 GeV Pb-|-Pb and 
V^Af7V~200 GeV Au-f An collisions, constrain the curva- 
ture parameter within a narrow range, k ~ 0.018 — 0.023. 
It is factor of ~ 3 larger than in the lattice QCD calcu- 
lations. To our knowledge, this is the first attempt to 
obtain QCD critical line from experiment data. 

II. MODEL FOR J/i/> SUPPRESSION IN QGP 

In a deconfined medium (QGP), inter-quark (color) 
potential is screened. The screening increases with tem- 
perature and above a certain temperature, the poten- 
tial is too weak to bind charm and anti-charm quarks 
into a bound state. Charmonium states are then sup- 
pressed. Matsui and Satz [l^, first predicted the phe- 
nomena. Over the years, in several experiments, J/ip 
yield in heavy ion collisions has been measured. Measure- 
ments at SPS (yi^^=17.3 GeV) and RHIC (v^^^=200 
GeV) energy validate Matsui and Satz's predictions. 

To mimic the onset of deconfining phase transition 
above a critical energy density and subsequent melting 
of JjA/j'SjOGP motivated threshold model was proposed 
in [1^, [2^. In the threshold model, Jf-i/j suppression 
is linked with the local transverse density. If the local 
transverse density, at the point where J/tjj is formed, 
exceed a critical or threshold value, J/ip^s are melted. 
In the present paper, we use a variant of the threshold 
model, if the local temperature of the fluid, exceed a 
threshold value Tj/^, J/t/^'s are melted. Approximately 
30-40% of observed J/^'s are from feed-down decay of 
the higher states, -0' and x [2lJl. Higher states, ip' and 
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FIG. 1: (color online) Lattice QCD calculation of QCD phase 
diagram, in imaginary chemical potential method and Tay- 
lor expansion method [1^1 a,re shown. The red line is the 
chemical freeze-out curve obtained in a statistical model [l7j . 
The phase diagram obtained in the present analysis is shown 
as the blue circles. 



X are expected to melt at temperature T^',T^ < Tj/^. 
Lattice QCD calculations indicate that J/ip^s can sur- 
vive up to Tj/^ — 1.5 — 2Tc [2l|, higher states survive 
up to T^i « « 1 — 1.2Tc. For sequential melting 
of charmonium states, if the melting temperatures Tj/^ 
and T^',T^ are substantially different, J/iIj survival prob- 
ability should show at least two step distribution (see for 
example Fig. 4 in [32]). No such two step distribution is 
seen either in y/sj^^— 200 GeV Au-t-Au or in -^5^^=17. 3 
GeV Pb-|-Pb collisions. Absence of two step structure in 
J ftp survival probability indicate that the melting tem- 
peratures of J/ijj^s and higher states are not substan- 
tially different. In the following we assume direct and 
feed-down J/i/j's are melted at the common temperature 
= KTc. In a baryon rich medium, melting temper- 
atures Tj/^ depend on the QGP phase diagram through 
the critical temperature T'c(/x) as in EqlT] 

Threshold model require the local temperature (Ti) 
and (baryonic) chemical potential (/i^) of the fluid pro- 
duced in nuclear collisions. We propose to obtain them 
from experimentally measured quantities, e.g. charged 
particles rapidity density ^^y^ and net baryon rapid- 
ity density (B ~ B)^. Experimentally measured values 



of 



.0 m and {B~B)my=o [23l,[2 



0-5% 



Y^^^=17.3 and 200 GeV collisions, are listed in tableHl 
The initial entropy density is then determined from the 
well known relation [23, 



1.5 X -i.QdNch, 
T:R?Ti drj 



(2) 



In Eql21 Ti is the initial time, for which we have used 
the canonical value Ti=1.0 fm. We do note that Eql2]is 
obtained under the assumption of isentropic, longitudinal 



TABLE I: Rapidity density 



net baryon density {B — 
^jvjv^^'-'-' iiJiiL. ciiiu Y'^]vjv^200 GeV 
Au-|-Au collisions. Corresponding average temperature {Ti) 
and ^li are also listed. 



B)^ in 0-5% Vsjvjv = 17.3 Pb-^Pb and 



(GcV) 




(B-B)^|,=o 


(MeV) 


(MeV) 


17.3 (Pb+Pb) 


308 


75 


191.5 


390.2 


200.0 (Au+Au) 


690 


7 


242.4 


23.6 



expansion. Similarly, initial baryon density is obtained 
as. 



[B 



dn 



j/=0 



(3) 



(..) in Eql^JEl indicate that the entropy density and 
baryon density in Eas l2l3l must be understood as spa- 
tially averaged over entropy /baryon density. To intro- 
duce spatial dependence, we assume, that in an impact 
parameter h collision, initial entropy density and baryon 
density in the transverse plane is distributed as (26j , 



Si{x,y) = sq 
ni{x,y) = no 



(l-/) ^"'"'fr^^ +/^con(x,y) 



,(4) 
(5) 



So and no is fixed to reproduce the average entropy den- 
sity {s)i and baryon density {n)i, as determined from ex- 
perimental values. Npart{x,y) and Ncou{x,y) in Eq |4l5l 
are the transverse profile of the participant density and 
binary collision number, obtained in a Glauber model. / 
is the fraction of hard scattering. It is approximately con- 
stant in the energy range -y/s^jy^ 17-200 GeV, / sa0.12. 

Local entropy density {si(x,y)) and baryon density 
(ni{x,y)) can be converted into local temperature and 
baryonic chemical potential using an equation of state. 
In [27], a parametric form of QGP equation of state was 
given, which matches with the lattice QCD simulation at 
zero chemical potential and to the known properties of 
nuclear matter at zero temperature. Entropy density (s) 
and baryon density (n) in parameterised form are. 



4^ 
90 

Nf 



16 



21V/ 



Nf , 



2CT 



81tt' 



(6) 
(7) 



Nf w 2.5, 2C w 0.24, D w 0. In tableHl the average 
initial temperature and baryonic chemical potential in 
^^^=17.3GeV Pb-^Pb and ys^jv=200 GeV Au-hAu 
collisions are listed. From SPS to RHIC, chemical poten- 
tial decreases by a factor of ^15. Temperature however 
is increased by ~25%. 
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FIG. 2: (color online)The filled circles are PHENIX data [IH 
for nuclear modification factor in y^ijviv=200 GeV Au+Au 
collisions. Black lines (from bottom to top) are j/ip survival 
probability in the present model for k=0.03 and J/ip melting 
temperature Tj/^ = KTc, K=1.5, 1.55, 1.6, 1.65 and 1.7, 
Tc=0.16. The dashed blue lines are prediction obtained with 
fixed Tj/^ = 1.65rc and k=0 and 0.1 and. The dotted line is 
the suppression due to the CNM efTect. 



To obtain J/ip survival probability {Sj/^), we ran- 
domly distribute J/ip^s in the reaction plane and count 
the number of J/tp^s that survive the threshold condition, 
Tiocai < Tj/^. We note that apart from suppression in 
deconfined matter, J/ip's are suppressed in cold nuclear 
matter also (CNM effect). It is important to eliminate 
the suppression due to CNM effect. If cold nuclear mat- 
ter effects are accounted for, the survival probability of 
J/ip's can be obtained as, 

Sj/ip = SqgP X ScNM (8) 

where Scnm is the survival probability in cold nuclear 
matter and Sqgp is the survival probability in QGP. 
In ^/sjY^=200 GeV d-|-Au collisions, J/tp production is 
consistent with cold nuclear matter effect quantified in a 
Glauber model of nuclear absorption with aats = 2 ± 1 
mb (Hi. aabs « 7.6 ± 0.07 mb ^22,] in 75^^^=17.3 GeV 
Pb+Pb collisions. 



III. RESULTS 

The present model has two parameter, the threshold 
temperature for J/tp melting, Tj/^ = KTc and the cur- 
vature parameter k of the QCD critical line. The crit- 
ical temperature is fixed to Tc(/i = 0)=160 MeV. The 
free parameters are obtained by fitting experimental data 
on J/tp production in -^5^^^=17.3 GeV Pb-|-Pb (29i] and 
yi^^=200 GeV Au-hAu [Ml collisions. In y/sj^j^=200 
GeV Au-|-Au collisions, the fluid is essentially baryon 
free and model predictions do not depend sensitively on 



FIG. 3: The solid circles are NA60 data for measured J/tp 
over expected J/tp in 158 AGeV Pb-|-Pb collisions. The solid 
lines (top to bottom) are present model predictions with 
K=0.014, 0.016, 0.018, 0.020, 0.022 and 0.024 respectively. 
Tj/^ = l.eOTc, Tc=160 MeV. 



the curvature parameter. The dependence on melting 
temperature however is strong. In FiglJl for thresh- 
old temperature Tjf^/Tc=1.5, 1.55,1.6, 1.65, and 1.7, 
model predictions are compared with the PHENIX data 
[3l| for nuclear modification factor in y^jy^=200 GeV 
Au-|-Au collisions. CNM effect is quantified in Glauber 
model of nuclear absorption with CTafcs^S mb. Exper- 
imental data are well explained for the threshold tem- 
perature ^^r^ = 1.60 ± 0.03. The melting temperature 
Tj/,i, « 1.60rc is consistent with lattice QCD calculation 
|21j . The blue lines in Figl2] are model predictions with 
K=0.0 and 0.1, Tj/^, = 1.6Tc. They cannot be distin- 
guished. 

Curvature parameter is constrained when J/tp^s are 
absorbed in baryon rich plasma. In 158 AGeV Pb-|-Pb 
collisions fluid is baryon rich and J/ip suppression de- 
pend sensitively on the curvature parameter k. In FiglSl 
centrality dependence of J/tp suppression in ^/spfj^— 17.3 
GeV Pb-I-Pb collisions is shown [29], [30]. The data are 
shown as the ratio of measured J/tp^s over the expected 
J/4>'s (in a Glauber model of nuclear absorption). Upto 
Npart—'^OQ, the ratio is consistent with unity. Suppres- 
sion in addition to the CNM effect is required in Npart > 
200 collisions only. In Figl3l the solid lines are thresh- 
old model predictions with Tj/^ — 1.60Tc and k=0.014, 
0.016, 0.018, 0.020, 0.022 and 0.024 respectively. Unlike 
in -y/s^j^=200 GeV Au-|-Au collisions, in -^5^^^=17.3 
GeV Pb-|-Pb collisions, J/tp survival probability shows 
sensitive dependence on k. For k < 0.016, J/tp^s are 
not suppressed in QGP medium. Best description to the 
data is obtained with k «0.02 {x^/N « 1.0). Description 
to the data deteriorate for lower or higher curvature. If 
uncertainty in the melting temperature is accounted for, 
Pb-|-Pb data fix the curvature parameter within a narrow 
range, k = 0.018 - 0.023, for Tj/^ = (1.57 - 1.63)Tc. 
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In FiglH QCD critical line obtained with the curvature 
parameter k — 0.018 — 0.023 is shown (the blue circles 
with error bars). Curvature k = 0.018 — 0.023, obtained 
from J/^ data is ~ 3 larger than in lattice QCD sim- 
ulations (k ~ 0.006). It is also closer to the chemical 
freeze-out curve. The result is encouraging with respect 
to detection of QCD critical end point. If QCD criti- 
cal line is close to the chemical freeze-out, possibility of 
distortion of signals associated with CEP, e.g. charge 
fluctuations etc., will be lessened. 

Before we summarise the results some limitations of 
the model are noted. Experimental information of parti- 
cle density and baryon density is converted to local tem- 
perature and baryonic chemical potential using the EOS 
(see Eq l6|7p . The EOS is for massless quarks and glu- 
ons, adjusted to reproduce zero chemical potential lattice 
QCD simulations. How far it represent the real world is 
uncertain. The model is also a static model. It neglect 
the possibility that J/ip's, initially produced in a cold re- 
gion and survived, at a later time can drift into a hotter 
region and get melted. Such possibilities, if accounted for 



would reduce the survival probability. Consequently, the 
present model give a lower limit of curvature parameter 

K. 



IV. SUMMARY AND CONCLUSIONS 

To summarise, we have shown that in a threshold 
model, J/ip suppression is sensitive to the curvature pa- 
rameter of the QCD phase diagram. In threshold model, 
J/ijj suppression depend on the local temperature and 
chemical potential. A model is proposed to obtain local 
temperature and chemical potential, from experimentally 
determined quantities e.g. rapidity density, net baryon 
rapidity density. QCD critical line with curvature pa- 
rameter K=0. 018-0. 023 is consistent with experimental 
centrality dependence of J/ip suppression in y^^^=17.3 
GeV Pb+Pb and ^/sj^j^=200 GeV Au-^Au colhsions. 
K=0. 018-0. 023 is approximately 3 times larger than that 
obtained in lattice QCD simulations. 
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